Theory of reduced superfluid density in underdoped cuprate superconductors 
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The critical temperature of an underdoped cuprate superconductor is limited by its phase stiffness 
p. In this article we argue that the dependence of p on doping x should be understood as a 
consequence of deleterious competition with antiferromagnetism at large electron densities, rather 
than as evidence for pairing of holes in the a; = Mott insulator state, p is suppressed at small 
X because the correlation energy of a d-wave superconductor has a significant pairing-wavevector 
dependence when antiferromagnetic fluctuations are strong. 
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I. INTRODUCTION 



phase stiffness of a superconductor: 
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The fascinating and rich phenomenology of high tem- 
perature cuprate superconductors has been very thor- 
oughly studied over the past 20 years. Although there 
is substantial variability in detail from material to ma- 
terial, all cuprates exhibit robust Mott insulator anti- 
ferromagnetism when the hole-doping fraction x is very 
small, superconductivity which appears when x exceeds 
a minimum value ~ 0.1, and a maximum T^, in optimally 
doped materials with x 0.2. In the underdoped regime, 
the superconducting transition temperature is limited by 
phase fluctuationai'^'^'^, and experiments hint at a wide 
variety of (typically) short-range correlations associated 
with competing charge and spin orders. The underdoped 
regime poses a fundamental challenge to theory because 
its electronic properties are not fully consistent with any 
of the various well-understood fixed-point behaviors that 
often help us to classify and predict the properties of very 
complex materials. 

The phenomenological parameter p used to character- 
ize phase-fluctuation stiffness in a superconductor is nor- 
mally expressed in terms of the superfluid density Ug by 
writing p = h^ris/m* , an identification that is partly jus- 
tified by BCS mean-field theory. The increase of p with x 
in cuprate superconductors is therefore readily accounted 
for by theories"' in which superconductivity is due to 
the condensation of Cooper pairs formed from holes in 
a doped Mott insulator—. Theories which start with this 
view must still explain the fact that p vanishes at a non- 
zero value of x, and deal with the awkward property that 
cuprate superconductivity evolves smoothly from the un- 
derdoped regime to an overdoped regime in which it ap- 
pears to be explainable in terms of conventional band- 
quasiparticle Cooper pair condensation. In this article we 
propose an alternate explanation for the ^-dependence of 
p based on band-quasiparticlc pairing. Our argument ac- 
counts for the correlation energy of a d-wave supercon- 
ductor in the presence of incipient antiferromagnetism 
and is based on the following general expression for the 



P = 



1 d^E 



(1) 



where A is the area of the system, P the pairing 
wavevector^, and E is the total energy including both 
mean- field and correlation contributions: E — E^^^ 
j^cor^ The familiar BCS theory expression for p captures 
only the mean-field theory contribution to the energy. 

When superconductivity is viewed as a weak-coupling 
instability of a Fermi liquid, it is usually implicitly as- 
sumed that E'^"^ is not significantly infiuenced by the 
formation of the superconducting condensate, and cer- 
tainly not by changes in the condensate's pairing mo- 
mentum P. In the case of simple models with parabolic 
bands and galilean invariance, neglect of the correlation 
energy contribution can be justified rigorously. We ar- 
gue the correlation energy contribution is significant in 
underdoped cuprates because there is direct competition 
between the Fermi sea quantum fluctuations which con- 
dense in antiferromagnetic and d-wave superconducting 
states. Consequently the pair-breaking effects of finite 
P, which weaken superconductivity, also increase the im- 
portance of antiferromagnetic fiuctuations, lowering E'^°^ 
and decreasing p compared to its mean-field value. In 
the following sections we first use a fully phenomenolog- 
ical and then a partially microscopic extended-Hubbard- 
model weak-coupling theory to expand on this idea. The 
conjugate relationship- between pairing and antiferro- 
magnetic fluctuations plays an important role in the fluc- 
tuation spectrum and hence in theories of the correlation 
energy. In our theory of the underdoped state, the reso- 
nant magnetic mode (INSR) observed in inelastic neutron 
scattering-'-?^ experiments therefore has a somewhat dif- 
ferent interpretation than in most earlier theor y^^i^^i^^i^^ , 
appearing as a kind of magnetic plasmon. 
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II. PHENOMENOLOGICAL THEORY 

The basic ideas of our theory are quahtative, inde- 
pendent of most microscopic details, and most easily 
described in terms of the properties of a low-energy 
effective-field model for the collective fluctuations of a 
weak-coupling d-wave superconductor. The relationship 
to less transparent generalized random-phase approxima- 
tion (GRPA) correlation energy calculations is explained 
below. We construct a quantum action by introducing 
a set of states which incorporate the coupled triplet- 
pairing and spin-density fluctuations on which we focus. 

y]) is the Fock-space Slater determinant ground 
state of the quadratic Hamiltonian 
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(For notational simplicity we have exhibited here only 
fluctuations with zero spin projection along the quan- 
tization direction.) In Eq.[[2], t labels the four neigh- 
bours of each site on a two-dimensional square lattice, 
and {—Y represents the d-wave variation of mean- field 
near-neighbor pair potentials. Using these states as an 
approximate identity resolution leads to the following 
low-energy imaginary-time action for the collective vari- 
ables 6i and Vi: 



S = 



dT 



h {n<l>.v]\dAn<l>.v]) + E[^,V] , (3) 



where E[(I),V] = 1/]) and H is the fuU 

microscopic Hamiltonian. Mean-field theory states are 
obtained by minimizing E\(t>i The first term in the ac- 
tion captures the Berry phase coupling^ between pairing 
and spin-density fluctuations on which we now elaborate. 

The potentials associated with the two types of fluctu- 
ations are: 



The Berry phase term can be evaluated explicitly for 
small fluctuations by using perturbation theory expres- 
sions for the wavefunctions which appear in the Slater 
determinant | \1/ [(/), F ] ) : 



/•OO 

Jo - 



(5) 



where 



2E- 



Im 



(6) 



d-wave mean-field state of a generalized Hubbard model with 
U/t = 2.0, V/t = 2.0, t'/t = -0.3, and x = 0.12. A 
(ci|Cj|) = (— )'^0.145 in this case. 



In Eq.( [6]) we have made the usual Nambu spin-down 
particle-hole transformation of the mean-field Hamilto- 
nian so that it has two eigenstates at each wavevec- 
tor in the square lattice Brillouin zone with eigenval- 
ues iEp, one ixp,-) occupied and one (xp,+) unoccu- 
pied. In Fig.[[T] we show Berry curvature values calcu- 
lated from this expression as a function of k which are 
strongly peaked near k = Q — (7r/a,7r/a); these results 
arc robust over a broad range of dopings, gap sizes, and 
band-structure models. Pairing phase fluctuations are 
conjugate to spin-density fluctuations for k near Q, just 
as they are conjugate to charge-density fluctuations for 
k near 0, because of^ the d-wave property '^^^q — ~^p- 

We now argue that there is competition between the 
correlation energy gain due to antiferromagnetic fluctua- 
tions and d-wave singlet Cooper pair formation. Strong 
experimental evidence for this competition is provided by 
the apparent enhancemen't^'i^ of antiferromagnetic fluc- 
tuations in cuprate vortex cores and in cuprates placed 
in an external magnetic fields. Changes in external con- 
ditions which weaken superconductivity enhance antifer- 
romagnetism. Here we explore consequences of this com- 
petition for the correlation contribution to the superfluid 
density, i.e. for the dependence of correlation energy on 
pairing-momentum P. 

In our model the quadratic fluctuation action of a d- 
wave superconductor is 



(7) 



In Eq.[[7], K't' and K'''' are phase and spin-density stiff- 
nesses. The onset of antiferromagnetism occurs when 
— 0. In using this action we assume that the 
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most important quantum fluctuations are d-wave pair 
phase and spin- fluctuations. The microscopic calcula- 
tions GRPA described in the next section support this 
assumption. In Eq.[[7], frequency dependence is indi- 
cated in Crr, , and Ki to recognize the existence of 

^ k k 

non-adiabatic effects accounted for in these more micro- 
scopic calculations but neglected in this qualitative dis- 
cussion. The quadratic fluctuation action then describes 
a system with collective modes at energies 



k k 



Ck 



and a corresponding zero-point energy contribution 



(8) 



(9) 



This adiabatic theory of the INSR mode is accurate only 
when lies below the particle-hole continuum; the 

prime on the wavevector sum above recognizes that this 
condition is satisfied only near k = Q. The fluctuation 
correction to p can be related to the pairing-wavevector 
dependence of the zero-point energy as we explain below. 

We expect K''p{Q) to decrease with P because sup- 
pressed pairing favors antiferromagnetism. The strength 
of this dependence can be estimated roughly from ex- 
periment by associating the magnetic length £b at the 
magnetic field strength required to induce antiferromag- 
netism in a cuprate superconductor with the value of 
P~^ at which K^^IQ) goes to zero. Taking a typi- 
cal value for this field ^ 100 Tesla and assuming that 
the resonance mode is well defined over the portion of 
the BZ with large Berry curvature (say ~ 10%), gives 
a negative correlation energy contribution to the phase 
stiffness per two-dimensional cuprate layer of p'^"'^ ~ 

-0.1n£| E''^" ~ 0.05eV, comparable to the 

value of p inferred from penetration depth measurements 
in optimally doped samples. Although this estimate is 
clearly very rough, it does establish that the correla- 
tion correction to p can be substantial in the underdoped 
regime. 

Charge density fiuctuations are not included in this 
analysis because their Berry phase coupling to the phase 
fluctuations of the superconducting order parameter is 
large only near fc = and negligible near k = Q^. 
Thereofore the charge density fluctuations do not play 
a signiflcant roles for the physics near k = Q, instead 
they cause the instabilities near fc = in the microscopic 
GRPA calculations which we mention below. 



III. MICROSCOPIC GRPA THEORY 

We now evaluate the correlation energy of an extended 
Hubbard model^'' in the generalized random phase ap- 
proximation (GRPA) approximation. The model we 



study in this article has on-site repulsive interactions U 
which drive antiferromagnetism and near-neighbour at- 
tractive interactions V which drive d-wave superconduc- 
tivity: H = Ht+Hu +Hv, 



Ht^-t J2 clc.a + h.c. - 1' J2 



Cja + h.c. 



<i,j>,(y 



<i,j>' ,cr 



-V ^ fliafljo 

<i,j>cra' 
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In Eq.( \TUI V, t and t' should aU be thought of 
as effective parameters which apply at the energy scale 
of pairing and depend on x. Values for V, t, and t' 
can be estimated from ARPES data^. Spin dependent 
Heisenberg near-neighbor interactions of the type used in 
t — J models could also be used in the low-energy effec- 
tive Hamiltonian, but are neglected here for simplicity. 
The conclusions we draw in this article do not depend 
on whether the near-neighbor effective interaction which 
drives d-wave superconductivity is spin-independent or 
spin-dependent . 

The GRPA correlation energy of a d-wave condensate 
state with pairing momentum P is^^< 



ph 



(Ag)], (11) 



where e^''(P, q) is a quasiparticle particle- hole exictation 
energy, and uji {P, q) is the corresponding GRPA excita- 
tion energy. This equation can be derived by expanding 
the GRPA Hamiltonian to quadratic order in particle- 
hole excitation amplitudes approximated as independent 
bosons. The correlation energy expression then drops 
out of a boson Bogoliubov transformation. The analysis 
in the preceeding qualitative discussion assummed that 
the P-dependence of E'^°^ is dominated by its collective 
mode contribution, an assumption that is largely justifled 
by the following more microscopic calculation. 

In a GRPA theory excitation energies uJi [P, q) are 
obtained from a time-dependent mean-fleld-theory^^ in 
which the quasiparticles respond to the external potential 
and to induced mean-fleld potentials: H' — H^^*^ + (t) 
where Hf{t) = Hf\t) + Hf^{t), 

p,k.q,a 

+G{k,p, q) [(5(ct+g-„ •=P'^)4-^,-a ^ka 
+Hik,p) [S{cl_.^ cl^)cj:^ c._-,^ + h.c. 



p,k.q,o' 



+J{k,p) ['5(ct._p-^ct^)cj^Vfe^ + 

(12) 
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where 

F{q) = U — 2T^[cos + cos qy] , 

G{k,p, q) = 2V[cos{kx -Px - qx) + cos(fcy - Py - qy) 

- cos qx - cos qy] , 
H{k,p) = 11/2 - "l^[cos(fc2, -p.j,) + cos{ky - Py)], 
I{k,p) = [-[/ + 2V{cos{ka: ~Px) +cos{ky - Py))], 
J{k,p) = -y[cos(fcj; -p,^) +cos(fcy -py)], 

(13) 

and a = —a. Because the d-wave BCS ground state is 
a spin-singlet, its elementary excitations consist of = 
1 triplet and S = Q singlet branches. Hf^{t) captures 
the Sz = ±1 portions of the triplet fluctuations, studied 
by Demler and Zhang in a different context-S. Hf^{t) 
captures singlet and triplet Sz = fluctuations. 

Applying linear response theory to the quasiparticle 
Hamiltonian H — Hmf{P) + H'{t), we can compute the 
change in element of density matrix S{pab) by: 

5{Pab{t)) ^jj dt'9{t-t') {[H'{t'),p,b{mH.,,, (14) 

, leading to an equation of the form: 

{uji-M)p{q,uj) = H^-\ (15) 

where p is a column representing the change in the quasi- 
particle density matrix. The collective mode energies 
{ijji{q)} are the eigenvalues^^ of the matrix M which can 
be read off Eq. ( I14p . For a given pairing momentum P 
and excitation momentum (f, the number of particle hole 
pairs is proportional to the number of momenta in the 
Brillouin-zone, which we limit by using periodic bound- 
ary conditions with a finite quantization area A — L^. 
Diagonalizing M is equivalent to performing the boson 
Bogoliubov transformation, and equivalent to summing 
the ladder and bubble diagrams used to represent the 
GRPA in diagramatic perturbation theory. 

Typical GRPA extended Hubbard model results for 
the spin and pair response functions of the P — Q state 
at wavevector q = Q exhibit a single collective mode 
below the particle-hole continuum with large weight in 
both responses. In Fig. [[2], we illustrate the dependence 
on U of the energy of this excitation, and of its weight in 
spin-density and triplet-pair response functions, with V, t 
and t' (and hence the d-wave mean-field state) held fixed. 
As U increases, promoting antiferromagnetism, the col- 
lective mode excitation energy decreases, its weight in 
the spin response function increases, and its weight in 
the pair-excitation spectrum decreases. These proper- 
ties are consistent with our qualitative effective theory, 
given the expectation that K"-^ should decrease with U . 

The opposing variations of triplet-pair and spin-density 
weights demonstrates that the Berry phase mechanism 
dominates coupling between antiferromagnetic and pair- 
ing phase fluctuations as expected. 
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FIG. 2: Resonance mode energy and weights in spin and pair 
response functions vs. U for q = Q and doping x — 0.12. 
These results were obtained for P — and calculated with 
V/t = 2.0 and t' /t = -0.3 and 34 x 34 points in the BZ. For 
each value of U the vertical (horizontal) bar represents the 
collective mode weight in the spin (pair) response function. 
The weights were evaluated using the same response function 
defintions as Tchernyshyov et, al.^"^. 

IV. DISCUSSION 

Our qualitative discussion suggested that there should 
be a strong negative correlation contribution to p be- 
cause the INSR near q — Q softens with increasing 
\P\. Because we perform our calculations with peri- 
odic boundary conditions, we estimate the mean-fleld 
^pMF^ and correlation {p'^°^) contributions to p by com- 
paring energies at P = and P — P""" = (27r/L, 0); 
p « [£;(p™™) - E{Q)]/2tt^. Fig.[ E] illustrates the 
fluctuation-wavevector dependence of correlation contri- 
butions to p. As expected we find thafe^^ modes near 
q — Q soften, making a negative contribution to p. 
This dependence of collective mode energies on the pair- 
momentum of the superconducting condensate is unusual 
and is indicative of the microscopic competition between 
antiferromagnetism and d-wave superconductivity. This 
result contrasts strongly with the absence of any signif- 
icant dependence of plasmonic collective modes on pair 
condensate properties in conventional superconductors. 

As indicated in Fig.[ [3] we also find that for the 
model parameters chosen, collective modes at momenta 
near (0, 0) have complex energies for both values of 
P. This finding reflects the tendency of extended Hub- 
bard models, and of real cuprate materials, to longer 
period density-wave instabilities^. We do not believe 
that these ubiquitous instabilities, which appear to be 
material specific, should not play an essential role in 
underdoped-cuprate superfiuid density suppression since 
long-wavelength density-wave order will have little im- 
pact on near-neighbour antiferromagnetic fluctuations. 

We find that p'^°^ / p'^'^ is negative and of order —1 in 
the underdoped regime when extended Hubbard model 
parameters are in the range thought to represent un- 
derdoped cuprates. We conclude that a substantial 
suppression of the superfiuid density due to the pair- 
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ing wavevector dependence of the correlation energy oc- 
curs in underdoped cuprates and that it is responsi- 
ble for the downturn in the critical temperature. Our 
weak-coupling theory is unable to describe physics very 
near the termination of superconductivity on the under- 
doped side, although there is some indirect evidence from 
experimental^ (for example from the relatively weak tem- 
perature dependence of p) that critical fluctuations are 
important in a relatively narrow doping range. Our ex- 
planation for reduced superfluid density in underdoped 
cuprates is independent of the microscopic origin of the 
effective near-neighbor interaction responsible for V and 
hence for d-wave superconductivity. 



FIG. 3: (Color online) Correlation contribution to p from 
fluctuations with wavevector k {p'^°^{k) for the same model 
parameters as in Fig.[[T]. p'^°^{k) is normalized so that its 
Brillouin-zone average is p'^°''/p*^^, which has the value —1.6 
for these parameters. The GRPA excitation energies in the 
black area near k = (0, 0) are imaginary reflecting longer 
length scale instabilitiea^ of the extended Hubbard model 
we use. These long-wavelength instabilities are sensitive to 
model details and independent of the p'^"'" contributions from 
near A; = tt, vr. 
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